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ABSTRACT 

We present phase resolved optical spectroscopy and photometry of V4580 Sagit- 
tarii, the optical counterpart to the accretion powered millisecond pulsar 
SAX J1808. 4-3658, obtained during the 2008 September/October outburst. Doppler 
^L} | tomography of the N in A4640.64 Bowen blend emission line reveals a focused spot of 

, emission at a location consistent with the secondary star. The velocity of this emis- 

sion occurs at 324 ± 15 km s" 1 ; applying a "if-correction" , we find the velocity of 
the secondary star projected onto the line of sight to be 370 ± 40 km s _1 . Based on 
. existing pulse timing measurements, this constrains the mass ratio of the system to 

be 0.044l°:oo4i and tne mass function for the pulsar to be 0.44^^3 Mq. Combining 
this mass function with various inclination estimates from other authors, we find no 
evidence to suggest that the neutron star in SAX J1808.4— 3658 is more massive than 
k> the canonical value of 1.4 Mq. Our optical light curves exhibit a possible superhump 

modulation, expected for a system with such a low mass ratio. The equivalent width of 
■ the Ca II H and K interstellar absorption lines suggest that the distance to the source 

is ~2.5 kpc. This is consistent with previous distance estimates based on type-I X- 
ray bursts which assume cosmic abundances of hydrogen, but lower than more recent 
estimates which assume helium-rich bursts. 

Key words: accretion, accretion discs - binaries: close - pulsars: individual: 
SAX J1808.4— 3658 - stars: individual: V4580 Sagittarii - stars: neutron - X-rays: 
binaries 



1 INTRODUCTION 

Low mass X-ray binaries (LMXBs) are systems containing a 
low mass secondary (M2 < 1 Mq) and a compact primary, 
either a neutron star (NS) or a black hole (BH). These may 
be divided into two subclasses based on whether they are 
persistent or transient sources. The transient systems or X- 
ray Novae (XRNe) are typically characterised by short pe- 
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riods of heightened luminosity separated by long periods of 
quiescence. During an XRN outburst the X-ray luminosity 
increases by 10 4 - 10 6 , reaching a sizable fraction of the 
Eddington luminosity. The X-ray outburst is accompanied 
by an outburst at UV/optical/IR wavelengths dominated by 
reprocessing o f X-rays in the accret ion disc (for example, see 
the review bv lCharles fc Coell2006l ). 

Until recently, the vast majority of known XRNe were 
black hole systems. As a result, systematic comparison be- 
tween the accretion discs in BH and NS transients was 
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difficult. At the same time, there was little direct evi- 
dence for the link between accreting LMXBs and isolated 
millisecond pulsars that had been anticipated prev iously 
|Alpar et al.ll 19821 ; iRadhakrishnan fc Srinivasanllf 984) . The 
discovery of a new class of transient, accretion-powered mil- 
lisecond X-ray pulsars (AMSPs), simultaneously addressed 
both of these issues. They provide dramatic confirmation of 
the li nk between acc reting LMXBs and millisecond pulsars 
(e.g. IWiinandsll2006l ): indeed, by combining these observa- 
tions with measurements of burst oscillations and pulsations 
from other LMXBs, NS spin periods are now known for 20 
LMXBs. The reason why the pulsations are observed in this 
type of system, with a short orbital period, may be related 
to the low accretion rate (expected for such short periods), 
permitting accretion to the polar caps of the NS even in the 
presence of a relatively weak magnetic field. 

SAX J1808.4-3658 is the proto typical AMSP. Initiall y 
detected by the BeppoSax mission (|in 't Zand et al.1ll998h . 
subsequent Rossi X-ray Timing Explorer (RXTE) obser- 
vations detected coherent millisecond pulsations at a fre- 
quency of 401 Hz from a NS in a ~2 hour binary or- 
bit. It was immediately realised that this was the long 
awaited missing link in t he evolution of LMXBs to iso- 
lated millisecond pulsars (|Wiinands fc van der" Khsl Il998l ; 
IChakrabartv fc Morgan! fl99Sl ). Subsequent obs ervations re- 
veale d the optical counterp art at R ~16.1 mag l|Roche et al.l 
1 19981 ). [H omer et ail l|200ll ) observed the optical counterpart 
in the quiescent state and found it to be at R ~ 21 mag. 
Using recent observ ations from the Gemini South telescope, 
iDelove et all (|2008T ) measure average quiescent magnitudes 
of g ~ 21.7 mag and % ~ 20.6 mag. The 0.5 - 10 keV X-ray 
luminosity (L x ) of the system is ~2x 10 36 erg s _1 in outburst 
din 't Zand et al.ll 19981, 12OOIJ) and 5 x 10 31 erg s _1 in quies- 
cenc e Jcampana et al.|l2002f) , for an assumed distanc e of 2.5 
kpc l|in 't Zand et aljboOlf ). iGallowav fc Cummind (|2006l ) 
have calculated that the distance to the source is ~3.5 kpc, 
assuming that the observed type-I X-ray bursts are helium- 
rich. 

The quiescent optical light curves exhibit sinusoidal 
variability, the phasing of which indicates that the optical 
maxim um occurs when the secondary is dire ctly behind the 
pulsar i|Homer et al.ll200ll ; lDelove et al.ll2008l ). The observed 
optical flux is in excess of that expected assuming repro- 
cessing of the measured X-ray flux. It was proposed that 
the irradia ting flux was provided by the spin-down fl ux from 
the pulsar l|Burderi et all2~003l ; ICampana et al.ll2004T ). which 
was irradiating the face of the companion star, hence pro- 
ducing the observed modulation. This effect is similar to that 
obse rved in the binary m illisecond pulsar PSR B 1957+20 
(e.g. iRevnolds et~al]|2007h . 

To account for their high spin frequencies, the neutron 
stars in AMSPs are expected to have accreted a few tenths 
of a solar mass during the lifetime of the binary, and hence 
should be systematica lly more massive than the ca noni- 
cal value of ~1.4 M dThorsett fc Chakrabar"tvlll999l) . Re- 
cent quiescent observations ( Heinke et al. 20081 ) have hinted 
that the NS in SAX J1808.4-3658 may be massive (Mi > 
1.4 Mq). Confirmation that the NS is indeed massive could 
have important implications for efforts to constrain the NS 
equation of state. 

For AMSPs, the orbital period (P or b) and projected 
semi-major axis of the primary orbit (and hence the pro- 



jected primary velocity, K\) are accurately known from anal- 
ysis of the X-ray pulse arrival times. Therefore for these sys- 
tems, only the projected secondary velocity [K-i) and the bi- 
nary inclination angle (i) are required for the primary mass 
to be fully constrained. Unfortunately, measurement of the 
radial velocity of the secondary star in SAX J1808.4— 3658 
in quiescence is not feasible, because of the faintness of the 
secondary star. However, as lSteeghs fc Casaresl (|2002h have 
shown in the case of the LMXB Scorpius X-l, it is possi- 
ble to identify narrow Bow en blend emission line features 
(e.g iMcClintock et al.1 1 19751 ) formed in the secondary star 
by fluorescence of UV photons produced in the inner accre- 
tion disc, and use these to determine its radial velocity - 
even when the binary is X-ray bright. This is likely the only 
means by which we will ever be able to measure the mass 
of the neutron stars in AMSPs, and confirm their expected 
massive nature. 

In addition, monitoring of these AMSPs provides us 
with an excellent opportunity to study accretion discs dur- 
ing outburst and to compare the observed structure with 
quantitative models. For example, in the short period, small 
mass ratio system s, the disc is expected to be ellipti- 
cal and processing (|Whitehurstlll988l ; iHaswell et al.ll200ll ). 
Doppler tomography of these systems may provide evi- 
dence for such an elliptical disc. In addition, strongly X- 
ray irradiated dis cs are also susceptible to warping out of 
the o rbital plane fPringlc 1996|; iFoulkes. Haswell fc Murray! 
2006) , which shou l d also manifest itself spectroscopically. 
See lElebert et aLI (2008) for the first such study of the 
AMSP HETE J1900.1-2455. 

In this paper, we present phased resolved optical spec- 
troscopy and photometry of V4580 Sagittarii, the opti- 
cal counterpart to the accreting X-ray millisecond pul- 
sar SAX J1808.4-3658, obt ained during the 2008 outburst 
l|Markwardt fc Swank1l2008l ). The paper is organised as fol- 
lows: In §2, we describe the spectroscopic and photometric 
observations. An analysis of these data is presented in §3, in 
which we show how Doppler tomography of the Bowen blend 
N in A4640.64 line allows us to measure the mass function 
of SAX J1808.4— 3658. Finally, we discuss these results in 
§4, and present our conclusions in §5. 



2 DATA 

Our data consist of phase resolved spectroscopy and pho- 
tometry of SAX J1808.4-3658, obtained during the 2008 
September/October outburst. These observations were mo- 
tivated by the need to obtain spectroscopy in the region of 
He 11 A4686 and the Bowen blend near A4640 in an effort 
to measure the radial velocity of the low mass seco ndary 
star (Af 2 ~0.05 M g ; iBildsten fc Chakrabartvl I2OO1I ) . and 
hence place constraints on the mass of the pulsar. Although 
SAX J1808.4-3658 is the brightest of all AMSPs in out- 
burst, no optical spectra of the system in this state have 
been published to date. 

Fig. [1] shows the RXTE All-Sky Monitor (ASM) light 
curve of the outburst, with the dates of our observations 
marked. The observations at CTIO 1.0-m and 0.9-m tele- 
scopes on mjd 54748 were taken at the same time, but have 
been offset slightly here for clarity. 
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Figure 1. RXTE ASM light curve, with the dates of our optical 
observations marked. 



2.1 Spectroscopy 

A total of 39 optical spectra of SAX J1808.4-3658 
were acquired with the FOcal Reducer and low disper- 
sion Spectrograph (FO RS) 1 and FORS2 spectrographs 
l|Appenzeller et al"1 ll998), mounted at the Cassegrain foci of 
the 8.2-m Kueyen (UT2) and Antu (UT1) Very Large Tele- 
scopes, respectively, at Cerro Paranal, Chile. The detector 
on FORS1 consists of a mosaic of two 2000 x 4000 pixel 
E2V CCDs (with 15 um pixels) while FORS2 has a mosaic 
of two 2000 x 4000 pixel MIT CCDs (also with 15 [im pix- 
els). The output data were binned by a factor of 2 in both 
spatial and dispersion directions. 

On 2008 September 27 UT, 16 spectra were obtained us- 
ing the 1200B grism on FORS1 with a 0.7 arcsec slit (seeing 
0.8 - 1.1 arcsec) and 360 s exposure times, while a further 23 
spectra were obtained using the 600B grism on FORS2 with 
a 1 arcsec slit (seeing 0.6 - 1.0 arcsec) and 240 s exposure 
times. In all observations, the SAX J1808.4— 3658 spectrum 
was on CHIP1 of the FORS1 and FORS2 mosaics. Appro- 
priate bias, flat and comparison HgCd arc lamp frames were 
also taken. 

The data reduction process was identical for both the 
FORS1 and FORS2 data. The frames were firstly processed 
using the iRAtQ CCDPROC routines to remove instrumen- 
tal effects. The spectra were then optimally extracted using 
tasks in the iraf kpnoslit package. Wavelength solutions 
for FORS1 (FORS2) were found by fitting 2nd order cubic 
splines to 16 (23) lines in the arc spectra, giving RMS errors 
of —0.016 (0.036) A, and these solutions were then applied 
to the science frames. The resulting wavelength range was 
AA3620 - 5050 (3330 - 6350), with a dispersion of -0.70 



IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



(1.49) A pixel 1 and a resolution (measured from arc lamp 
lines) of -1.6 - 1.9 (4.3 - 4.8) A. 

Because only one arc lamp exposure was taken during 
the daytime for each dataset, night skylines were examined 
in the SAX J1808.4— 3658 frames to correct for instrumental 
flexure. For both the FORS1 and FORS2 datasets, a sim- 
ilar correction routine was employed. Firstly, the relative 
offsets between the wavelength scales for each of the indi- 
vidual spectra were found by cross-correlating each of the 
sky spectra against the last spectrum in the s et (preliminary 
temp late), using the RVSAO package in iraf |Kurtz fe Minkl 
1998). These offsets were used to shift each of the individual 
spectra onto the wavelength scale of the preliminary tem- 
plate, and the resulting spectra were averaged to produce a 
master sky template. The positions of the strongest skylines 
were examined in this master template, in order to com- 
pute an average absolute shift. The master template was 
shifted by this amount, so that all lines appeared at their 
correct wavelengths, and the final shifts required for each 
spectrum were found by cross-correlating the spectra against 
the shifted master template. These shifts were then applied 
to both the data and error bands of each spectrum. The re- 
quired shifts for the FORS1 data were within ±0.07 A for 
each spectrum. For the FORS2 data, the required shifts were 
larger, —0.5 - 0.7 A. The reduction procedure was repeated 
and confirmed using the FORS pipeline softwar^f] 

The data were then imported into the MOLLY spec- 
tral analysis package, where they were shifted to the he- 
liocentric rest frame, and re-binned onto a common veloc- 
ity scale of 49 km s _1 pixel -1 for the FORS1 data and 
95 km s" 1 pixel" 1 for the FORS2 data. The continuum was 
removed by dividing by a cubic spline function (after mask- 
ing the emission and absorption features) and subtracting 
unity. The normalized averaged FORS1 and FORS2 spectra 
are shown in Fig. [2] The o rbital phase was det e rmine d us- 
ing the X-ray ephemeris of I Jain. Dutta fc Paul (|2008h . We 
transformed the To (epoch when the pulsar is at the as- 
cending node) of the ephemeris from tdb to UT, using the 
xTime toojfl assuming that for the purposes of this work, 
tdb is equal to tt and also that the heliocentre and barycen- 
tre are identical. We also shifted To forward by one quarter 
of an orbital period, so that phase zero represents supe- 
rior conjunction of the pulsar i.e. when the pulsar is at 90° 
orbital longitude from the ascending node. The ephemeris 
we used was therefore To = 2450915.398669 hjd ut and 
Porb = 0.0839022785 d. 

Because the primary lines of interest were the Bowen 
blend and He n A4686, the data near these wavelengths 
were normalized separately, to ensure optimum continuum 
removal. 



2.2 Photometry 

In an effort to constrain the binary inclination angle via fit- 
ting t he light curve with a physical model of the light sources 
(e.g. iDelove et aT]|200Sh . we obtained phase resolved pho- 
tometry of SAX J1808.4— 3658 from several telescopes dur- 
ing the 2008 September/October outburst. On 2008 October 



2 http:/ /www. eso.org/sci/data-processing/software/pipelines 

3 http:/ /heasarc. gsfc.nasa.gov/cgi-bin/Tools/xTime/xTime. pi 
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Figure 2. FORS1 (bottom), FORS2 (middle) and combined (top) averaged normalized spectra obtained on 2008 September 27. The 
spectra have been shifted to the rest frame of the system, using the systemic velocity of —83 km s — 1 found in §3, and offset vertically 
for clarity. 



1 ut, we obtained 20 Gunn/SDSS i-band images using the 
2-m Faulkes South Telescope, at Siding Spring Observatory, 
Australia. The EA02 camera was used, which has 2048 x 
2048 pixels binned 2x2 into effectively 1024 x 1024 pixels, 
each of which is 0.278 arcsec pixel -1 ; the total field of view 
is 4.7 x 4.7 arcmin. 

On 2008 October 7/8 UT we obtained 40 Gunn r-band 
images from the 1-m SMARTS /Yale Telescope, operated by 
the SMARTS consortium, at Cerro Tololo Inter-American 
Observatory (CTIO), Chile. A further 31 Gunn r-band im- 
ages were obtained from the same telescope on 2008 October 
8/9 UT. The CCD system used for the observations was built 
by NOAO and Fermi National Accelerator Laboratory as 
part of the development of the Dark Energy Camera project 



(see www.darkenergysurvey.org and iDePov et al.1l2008l . for 
more details). The detector in the system is a thick, fully 
depleted 2048 x 2048 pixel CCD with 15 um pixels, which 
corresponded to approximately 0.3 arcsec pixel -1 on the sky. 
The exposure time for all observations was 300s. The raw 
images were bias corrected, trimmed, and flat-fielded using 
the CCDPROC routines in iraf. 

On 2008 October 9 UT we obtained 9 Bessel V-band 
images using the 0.9-m SMARTS Telescope, also at CTIO. 
On October 12, we acquired 8 Bessel J?-band images. The 
SITe CCD has 2048 x 2046 pixels, with 0.4 arcsec pixel -1 . 
As the CCD was read out using 4 amplifiers, the frames 
were bias corrected and flat-fielded using tasks in the I RAF 
quadproc package. 
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Figure 3. Light curves of SAX J1808.4-3658, folded on the 
ephemeris defined in §2.1, and plotted twice for clarity, (a) i-band 
data obtained with the Faulkes South telescope on 2008 October 
1 UT. (b) r-band data obtained with the CTIO 1-m telescope on 
2008 October 7/8 UT. (c) Bessel .R-band data obtained with the 
0.9-m CTIO telescope on 2008 October 12 UT. All three light 
curves have the same magnitude scale. 



Photometry was performed using the daophot ()Stetsonl 

1 19871 ) point spread function fitting package in iraf. 
The magnitude of SAX J1808.4-3658 was found by 
comparison with the thre e star s listed in table 1 of 
iGreenhill. Giles fc Couturesl (2006). For the i- and r-band 
frames, we firstly transfor med the .R-band magnitudes to 
r-band using equation 4 of Ijordi. Grebel fc Ammonl (|2006l ) 
(for V — R ^ 0.93) and /-band to i-band using equation 8 of 
Ijordi et al.l (|2006r ). Fig. [3] shows three light curves, phased 
using the ephemeris defined in H2.ll and plotted twice for 
clarity. Two of our data sets were obtained contemporane- 
ously, and these are plotted in Fig. [4] 



3 RESULTS 
3.1 Photometry 

The main purpose of our photometry was to identify a mod- 
ulation that could be used to constrain the binary inclination 
a ngle, by fitting the data usi ng, for example, the elc code 
of|o rosz fc Hauschild"tl(|2000l ). 



_ (a) CTIO 0.9-m 2008 October 9 7-band 



(b) CTIO 1-m 2008 October 8/9 r-band 
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Figure 4. Light curves of SAX J1808.4-3658. (a) V-band data 
obtained with the 0.9-m CTIO telescope on 2008 October 9 UT 
(circles), (b) r-band data obtained with the CTIO 1-m telescope 
on 2008 October 8/9 UT (squares). Note that the r-band light 
curve has been shifted down by 0.1 mag for clarity. 



Out first light curve was obtained by Faulkes South just 
as SAX J1808.4— 3658 was beginning to decline from the out- 
burst. Although the minimum occurs near phase zero, the 
overall variatio n is more reminiscent of a superhump modu- 
lation (see e.g. lPatterson et al.ll2005 | ), or indeed a pre- eclipse 
hump from the hotspot (see e.g. iGrauer et al"1ll994l ). than 
a modulation due to the heated secondary star. Our sec- 
ond light curve was obtained when SAX J1808.4— 3658 had 
declined further from outburst maximum: this light curve 
exhibits a modulation closest to what might expected from 
a heated secondary (see Fig. [3}, although a superhump mod- 
ulation cannot be ruled out. If this modulation is due to a 
heated secondary, a reasonable fit to the light curve (us- 
ing elc) can be achieved for any of t he range of inclination 
angles found by iDelove et al.l ( 20081 ) . Our remaining light 
curves do not exhibit any obvious orbital modulation. Hence 
we conclude that these data cannot be used to place reliable 
constraints on the orbital inclination. Our light curves cover 
just over one orbital period each, so it also possible that the 
observed modulations are due to fluctuations in the mass 
accr etion rate. 

Ijordi et all (|2006h have very little data for V ~ R ~ 0, 
but assuming that their transformations are valid in this 
regime, we transformed the 2008 October 7/8 UT CTIO 1- 
m r-band magnitudes to Johnson- Cousins .R-band by using 
the simultaneous V-band data, to give V — R = 0.12 ± 
0. 04 mag. This is consist ent with the measurements in fig. 2 
of IGreenhill et al.l i|2006T ). De-reddening the R- and V-band 
magnitudes using our measured Av — 0.8 ± 0.3 (see H3.2|l 
yields (V - R)o = -0.08 ± 0.09 mag on this night. 
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Table 1. Equivalent width and full-width half-maximum for the 
strongest features in the SAX J1808.4— 3658 average spectrum. 



EW FWHM 

(A) (A) (km s- 1 ) 



He II A4686 


—0.85 ± 0.16 


23 ± 6 


1500 ± 400 


Bowen blend 


-0.35 ±0.16 






H/3 


5.0 ±0.4 


70 ±6 


4300 ± 400 


Ht 


2.8 ±0.3 


37 ±3 


2600 ± 200 


HS 


4.5 ±0.3 


48 ±3 


3500 ± 200 


He 


2.5 ±0.2 


32 ± 3 


2400 ± 200 


H< 


2.3 ±0.2 


29 ±3 


2200 ± 200 


H77 


1.2 ±0.2 


22 ±4 


1700 ± 300 


He I A4472 


1.0 ±0.3 


40 ± 15 


2700 ± 1000 




Interstellar feat 


ures 




Ca II K A3933.67 


0.77 ±0.07 






Ca II H A3968.47 


0.56 ± 0.08 






Na D a A5889.95 


0.97 ±0.12 






Na D 2 A5895.92 


0.68 ± 0.12 






DIB A5870 


0.14 ±0.04 







Nn/E (B -V) = 5.0 x 10 21 cm -2 mag" 1 jSavage fc Mathisl 
Il979 | . and referen ces therein), with a la error of 30% 
l|Bohlin et allll978l 1. this implies E(B-V) = 0.26^ og, con- 
sist ent with the value found u sing the DIB. 

unari fc Zwitterl (|l997h examine the relationship be- 
tween the EW of the Na Di line and E(B - V). The val- 
ues of E{B — V) we calculate above are only consistent 
with their analysis if the Na D absorption lines towards 
SAX J1808.4— 3658 have multiple components. Our spectral 
resolution is too low to determine if this is the case. The ra- 
tio of Na Di to Na D2 is ~1.4, which is also consistent with 
a relatively low E(B — V). 

Taking E(B - V) = 0.27 ± 0.11 A v = 0.8 ± 0.3. T his 
is consistent with the value found by I Wang et all l|200ll) . 

iMegier et all l|2005h established relationships between 
the EW of the Ca 11 K and H interstellar absorption lines 
and the distance to the source. Such a relationship suggests 
that Ca 11 is homo geneously distribute d throughout the in- 
terstellar medium l|Galazutdinovll2005l ). These relationships 
suggest a distance to SAXJ1808.4-3658 of 2250 ± 230 pc 
based on the Ca 11 K line and 2700±400 pc based on the Ca 11 
H line. The combined distance estimate is 2500 ± 400 pc. 



3.2 Spectroscopy 

Our spectra (Fig. [2} are dominated by broad Balmer ab- 
sorption lines, with interstellar Ca 11 H and K and Na D 
absorption lines, and emission lines from the Bowen blend 
near A4640 and from He 11 A4686. The diffuse interstellar 
band (DIB) at A5780 is also present. The main characteris- 
tics of these features are summarised in Table [1] 

The Balmer lines from SAX J1808.4-3658 itself extend 
to very high velocities, implying that they are produced in 
the optically thick accretion disc. There is some infilling of 
these lines by red-shifted emission features, most obviously 
in the H/3 absorption line. The H7 absorption line appears 
unusually narrow, compared with the other Balmer lines. 
It also appears to be blue-shifted relative to its rest wave- 
length, while the other Balmer lines seem to be relatively 
symmetrical about the rest wavelength, apart from the red- 
shifted emission. This may be due to higher levels of red- 
shifted emission in the H7 line, although why this should 
be is unclear. This effect is present in both the FORS1 and 
FORS2 data, and so would not appear to be a wavelength 
calibration issue. 



3.2.1 Extinction and distance 

E{B — V) is a function of the equivalent width (EW) of 
the DIB at A5870. Using four of the stars in table 1 of 
iJenniskens fc DesertJ || 19941 ) . EW(5780)/E(B — V) = 0.52± 
0.11. For SAX J1808.4— 3658 £W(5780) = 0.14±0.04, there- 
fore we fi nd E(B- V) = 0.27±0.11. This correlation is also 
shown bv lWebsterl (|l993T ) (fig. 5), and the point (0.27, 0.14) 
lies close to a best-fitting line through the low E(B — V) 
(< 0.5) points. The errors (here and in the remainder of 
the paper) represent the la uncertainties (or 68.3% confi- 
dence intervals for non-normal distributions) , calculated us- 
ing Monte Carlo simulations assuming Gaussian statistics. 

Fitting the X-ray spectra reveals that the hydrogen 
column density to the source (A^h) is equal to (1.3 ± 
0.1) x 10 21 cm" 2 jCampana. Stella fc Kenneall200Sl ). Taking 



3. 2. 2 Projected velocity of the donor star 

As discussed in SJT] for systems with low luminosity com- 
panions like SAX J1808.4— 3658, the only way to determine 
the projected velocity of the secondary star is by observing 
the radial velocity shifts of the Bowen emission lines near 
A4640. The primary features contained in the Bowen blend 
are N ill emission lines at AA4634.13, 4640.64, 4641.85 and 
4641.96, and C in emission lines at AA4647.4 and 4650.1. In 
the majority of systems so fa r observed, the N hi A 4640.64 
line is the strongest of these IjCornelisse et al.ll2008r ). 

In some systems, the velocity shift of these individ- 
ual lines can be observed as S-waves in a trailed spec tro- 
gram fe.g lSteeghs fc Casaresl [20021 ; ICasares et al.ll2006j ). In 
other cases, the S/N and resolution of the observations are 
such that the S-waves are not obvious, but the presence 
of velo city shifts is revealed by performing Doppler tomog- 
raph y ilMarsh fc Hornd Il988l; iMarshl l200lh of the spectra 
fe.g. ICasares et alj 120031 ) . ICasares et alj l|2004l ) report on a 
successful implementation of this technique for the AMSP 
XTE J1814— 338, which revealed the velocity of the d onor 
star emission to be 345 km s" - 1 . lElebert et aH i|2008l ) at- 
tempted this for another AMSP, HETE J1900. 1-2455, but 
the S/N of their data was insufficient to provide a reliable 
measurement. 

Before attempting to detect the secondary star in ei- 
ther the trailed spectra or Doppler tomograms, we firstly 
used the existing information about the system to deter- 
mine the expected range of projec ted secondary velocities. 
Analysis of the pulse a rrival times l|Chakrabartv fc Morgan! 
1 19981 ; Ijain et all 120081 ) provides us with very accurate val- 
ues of P orb (7249.156862 ±0.000005 s) and aisini (62.809 ± 
0.001 light-ms), which combine to give the primary velocity, 
Ki = 27raisini/Porb = 16.32056 ± 0.00026 km s"\ Since 
the primary is a NS, we assume an extreme primary mass 
range of 1.4 - 3.0 Mq, and that i lies between 5° and 90°. 

The mass function equation relates the mass of the pri- 
mary to these observable quantities. 
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where q is the binary mass ratio (= M2/M1 = K1/K2) By 
solving Equation [T] for the cubic in q, we calculated J\2 for 
a range of i and Mi, and plot 7\2 vs. Mi for several values 
of i in Fig. [S] This plot shows that, for Mi < 3 Mq and 
all possible inclination angles, K2 < 700 km s . Because 
the ephemeris is known so accurately from the pulse timing 
analysis, emission from the secondary should appear very 
close to the positive V y axis, and at a velocity less than 
700 km s~ , in our Doppler tomograms. 

In order to perform Doppler tomography, a parameter 
which must be known is the systemic velocity (7), as this 
determines the offset of any underlying S-wave from the rest 
wavelength. Ideally, this should be obtained using radial ve- 
locity studies of the secondary star: this is impossible here 
however. Instead, we determined a value for 7 by fitting a 
Gaussian to the He n A4686 line, after masking the line 
core. If the flux in the wings of this emission line originates 
in the inner accretion disc, then a fit to the wings of this line 
should provide a reasonable estimate of 7, where the flux is 
less contaminated by irregularities in the outer disc such as 
the gas stream/a ccretion disc impact region (although see 
the discussion by iMarshl 1 199Sl ) . For the NS LMXB Cygnus 
X-2 (|Elebert et al.ll2009l ). and BH LMXB XTE J1118+480 
(Elebert et al., in preparation.), the wings of this line, av- 
eraged over an entire orbital period, provide values for 7 
consistent (within the errors) with those found by analysis 
of the secondary star absorption line spectra. 

For the FORS1 data, the best fit gives 7 = -77 ± 
14 km s , while for the FORS2 data, we get a best fit 7 
of —77 ± 16 km s . Unfortunately, the Balmer absorption 
lines were too broad to confirm these values. 

Our Doppler tomograms were constructed using the 



maximum entropy method (MEM), as implemented in 
doppler. An initial value for 7 of —77 km s _1 was cho- 
sen for both the FORS1 and FORS2 tomograms. Since the 
N in A4640.64 emission line is generally the strongest in the 
Bowen blend, we first attempted to make Doppler tomo- 
grams of this line. The tomograms were iterated to minimize 
the reduced x 2 (xt) between the data and tomogram, while 
maintaining a smooth image. Both FORS1 and FORS2 to- 
mograms showed clear emission near the positive V y axis, 
where emission from the secondary is expected, and in both 
cases this emission was at a value of ~300 km s _1 . Using 
the 3<r range in systemic velocity, Doppler tomograms of 
the other Bowen blend lines did not show any emission from 
near the V y axis, nor any focused spots of emission at any 
phase or velocity. 

We then combined both datasets and a Gaussian fit 
to the wings of the averaged He 11 A4686 line gave a value 
for 7 of —83 ± 12 km s _1 . As an independent check on the 
value of 7, Doppler tomograms of the N ill A4640.64 emis- 
sion line were created for a range of 7 values, and the spot 
characteristics (peak and FWHM) were measured in each 
tomogram by fitting with a point spread function. When 
the ratio of the peak to FWHM was plotted against 7, a 
maximum was observed at 7 = —55 ± 20 km s~ , which we 
regard as consistent with our estimate from the He 11 A4686 
line. We note that the spot of emission is closest to the V y 
axis near 7 = — 80 km s - , and so we adopt a value for 7 of 
—83 ± 12 km s~ , as determined from the He 11 A4686 line 
fitting. 

We created several Doppler tomograms from these 
data using 7 = —83 km s _1 , with different velocity bin- 
nings, phase binnings and filtering. In all cases, we ob- 
served a bright spot of emission on the positive V y axis near 
300 km s _1 . Fig. [6] shows the Doppler tomogram for the 
case where we re-binnned all the data at 49 km s _1 , and 
into 20 phase bins. The centroid of the spot of emission is 
at a velocity of K eul = 324 ± 8 km s _1 . For Doppler tomo- 
grams within the la range of 7 values, the spot remained 
focused and close to the same value of 324 ± 8 km s" 1 . We 
estimate that the effect of the spot location dependency on 
7 increases the la error on K ern to 15 km s . Again, none 
of the other Bowen blend lines showed any focused spots of 
emission in their tomograms. The spectra extracted using 
the FORS pipeline software, and processed in an identical 
manner, yielded identical results, within the uncertainties. 

3.2.3 The K -correction 

ifem is the velocity of the inner hemisphere of the secondary, 
and a correction must be applied to obtain the velocity of 
the centre of mass of the secondary - the A"-correction. 
This correction is the ratio of the velocity of the Bowen 
emitting region to the velocity of the centre of mass of the 
secondary, and depends on how the secondary is irradiated 
by the cen tral source. By simulating the irradiation of the 
secondary, iMunoz-Darias. Casares fc Martmez-Paisl (|2005T l 
computed A'-correction values for a range of mass ratios, 
inclinations and disc flaring angles (the opening angle of the 
disc rim above the plane, a). They fit their simulation results 
with 4th order polynomials in q, and present the coefficients 
for various values of a, and for i — 40° and 90°. 

Since K\ is known to high accuracy from the pulse tim- 
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Figure 6. Dopplcr tomogram of the N in A4640.64 emission line 
using the combined FORS1 and FORS2 data. The Roche lobe 
of the secondary for K% = 370 km s — 1 and q = 0.044, the gas 
stream from the inner Lagrangian point and the velocity of the 
disc along the gas stream are overplotted. 



ing analysis, we estimated a value of q as the ratio K\ to 
K em , calculated a 7f-correction based on this, then iterated 
until a value of K2 was converged upon. We used the poly- 
nomial coefficients computed for i = 40° (although we found 
very little difference for the i = 90° coefficients), and cal- 
culated the if-correction for a range of values of a. The 
maximum value of a, for which irradiation of the secondary 
would still occur, was also calculated for the different val- 
ues of q. Based on this, a must be less than ~9°, since, 
regardless of the value of q, any value of a greater than 
9° would cause the secondary to be totally shielded from 
the primary/inner disc by the disc. Within the la range 
of values for K Bm (309 - 339 km s _1 ) and for all values of 
a < 9°, we find K 2 = 370 ±40 km s" 1 . Therefore, our 
estimate for K2 is quite conservative, as it takes the full 
range of possible values of a into account, since we have 
no way of determining a. Combining our measurement of 
K2 = 370 ± 40 km s~ x with the existing parameters and 
EquationEO we find q = 0.044^°^, f{M) = 0.44±g;i| M 
and Mi sin 3 i = 0.48t°;" Mq. 

3.2.4 He II X4-686 Doppler tomography 

Using 7 = —83 km s" 1 we created a Doppler tomogram 
of the He 11 A4686 emission line. The input data we used 
was similar to that used for the N ill A4640.64 line - ve- 
locity binned to 49 km s _1 pixel -1 , in 20 phase bins. Fig. 
[H] shows the Doppler tomogram of the He 11 A4686 line, 
with the Roche lobe of the secondary, gas stream velocity 
and velocity of the stream along the accretion plotted for 
K 2 = 370 km s" 1 and q = 0.044. The main feature of the 
tomogram is a bright spot of emission, along the gas stream 
trajectory. Whether this coincides with the location of the 




Velocity (km s 1 ) 



Figure 7. Trailed spectrum of the N III A4640.64 emission line 
using the combined FORS1 and FORS2 data, in 20 phase bins 
(top panel) and trailed spectrum computed from the Doppler to- 
mogram in Fig. [6] (bottom panel). 
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Figure 8. Doppler tomogram of the combined FORS1 and 
FORS2 data for the He II A4686 emission line. The Roche lobe of 
the secondary for K2 = 370 km s — 1 and q = 0.044, gas stream 
from the inner Lagrangian point and velocity of the disc along 
the gas stream are overplotted. 



gas stream/accretion disc impact point is unclear. There is 
also a roughly circular pattern, due to the accretion disc, 
although this is much more apparent in tomograms where 
the xt nas n °t been decreased to as low a value. Fig. 
shows the trailed spectrogram used to produce this Doppler 
tomogram, as well as the spectrogram computed from the 
Doppler tomogram. 



4 DISCUSSION 

4.1 The mass of the pulsar 

Fig. llOl shows a plot of K% vs. i for primary masses of 1.4, 1.8, 
2.2, 2.6 and 3.0 Mq. Based on our estimates for the system 
parameters, we calculate that for a primary mass between 
1.4 and 3.0 Mq, the inclination angle must lie between 33° 
and 44° for K 2 = 370 km s" 1 and between 29° and 51° 
for the lcr range of values in K 2 - For pulsar masses of 1.4 - 
3.0 Mq, our value of q leads to companion masses of ~0.06 
- 0.13 Mq. 

Our data provide us with a value for Mi sin 3 i — 
0. 48lo 14 Mq. In order to determine the primary mass, the 
only remaining parameter that must be found is the incli- 
nation angle. There have been several attempt s to constrain 
the or bital inclination of SAX J1808.4-3658. iDelove et all 
(2008) have published Gemini SDSS i- and g-band light 
curves of SAX J1808.4— 3658 in quiescence, and find that 
the 2cr range of values for the inclination is 32 - 74° , irre- 
spective of the primary mass. Hence, for an inclination of 
53° ± 11° (1(7), we derive a pulsar mass of l.Olg'4 Mq. 

X-ray measurements taken during the 2008 outburst 
also help to constrain the inclination angle. For example, us- 
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Figure 9. Trailed spectrum of the He II A4686 emission line us- 
ing the combined FORS1 and FORS2 data, in 20 phase bins (top 
panel) and trailed spectrum computed from the Doppler tomo- 
gram in Fig. 151 (bottom panel). 
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Figure 10. K2 vs. i for primary masses of 1.4, 1.8, 2.2, 2.6 and 
3.0 Mq. The dashed lines mark the allowed range of inclination 
angles if K2 = 370 km s — 1 , and the dotted lines show the range 
of angles at the 1<t limits on K2- 

ing observations from Suzaku, ICackett et al.l (|2009l ) fit the 
relativistically broadened Fe K line, and found a value of 
i = 55+1 de 

grees (90% confidence). Combining this esti- 
mate with our measurements leads to a primary mass of 
0.9 ± 0.3 M (Iff). By analysing RXTE d ata from the 
2002 outburst, llbragimov fe Poutanenl (|2009l) find a value 
of i = 60 ± 5°, consistent with the results of ICackett et al.l 
i|2009l ) and lDelove et al.l (|2008T ). within the uncertainties. 

Although all of these estimates are somewhat model 
dependent, it is clear that, when combined with our results, 
they do not support a heavy NS in SAX J1808.4-3658. This 
would suggest that the NS in SAX J1808.4-3658 has not 
accreted a significant amount of material, implying that the 
mass transfer has been non-conservative over the lifetime of 
the system. 

4.2 Accretion disc 

Our spectra show broad Balmer absorption lines, with some 
in-filling of the cores on the red sides of the lines. These ab- 
sorption lines originate in the optically thick accretion disc. 
The s pectra show simil arities to the spectra of KS Ursae Ma- 
joris (|Zhao et alj|2006h . an SU UMa type (i.e. exhibits su- 
perhumps during superoutbursts) cataclysmic variable sys- 
tem. In particular, the Balmer absorption lines show some 
infilling in both cases, with this effect more pronounced at 
longer wavelengths. The offset emission we observe may have 
an origin in an eccentric precessing accretion disc, as was 
also observed in th e spectra of the LMXB XTE J1118+480 
l|Torres et al.ll2003 ) ■ Helium absorption is at most very weak, 
which suggests that the accretion disc is dominated by hy- 
drogen. 

The optical photometry reveals a modulation which 
may be due to the phase dependent visibility of the irra- 
diated face of the secondary star, similar to the situation 



in quiescence (e.g. iDelove et al1l2008r ). However, it is more 
likely that this modulation is due to a superhumping accre- 
tion disc, or simply due to variations in the mass transfer 
rate onto the NS. Superhumps are periodic optical modu- 
lations, originally seen during superoutbursts of SU UMa 
dwarf novae (DNe), but a lso seen in some LMXBs in out- 
burst (jHaswell et al.ll200ll . and references therein). The su- 
perhump period (P s h) is typically a few percent longer than 
P or b, although we have insufficient data here to determine 
if this is the case for SAX J1808.4-3658. 

The Faulkes i-band light curve shown in Fig. EJa) is 
more like a sawtooth pattern, rather than the sinusoid-like 
variation expected if the modulation was due to the sec- 
ondary star. Although the minimum of this light curve oc- 
curs near phase zero, the maximum brightness is at phase 
0.7 - 0.8, rather than phase 0.5. The CTIO l-m r-band 
light curve in Fig. E^b) is more sinusoidal, but the am- 
plitude of the modulation is a factor of ~2 lower than 
in the i-band light curve. The CTIO 0.9-m i?-band light 
curve shows no discernible modulation. An other AMSP, 
HET E J1900.1-2455 exhibits superhumps (jElebert et all 
I2008T ). and in this case the amplitude of the superhumps 
varied by a factor of ~3 over 11 days, without any signifi- 
cant change in the X-ray brightness. For L x ~10 36 erg s _1 
the accretion disc is also likely to exhibit warping caused 
by the irradiation. Whether or not this is the case for 
SAX J1808.4— 3658 is unclear, but the combination of ap- 
sidal and warped precession may explain why some of our 
light curves appear to exhibit periodic variations, and others 
do not. We also note that most of our light curves cover a 
single orbital period only, and hence we cannot exclude the 
possibility that purely random variations in the mass trans- 
fer rate could give rise to the apparent superhump variability 
discussed here. 



5 CONCLUSIONS 

We have presented Doppler tomography of the N ill 
A4640.64 emission line for SAX J1808.4-3658, from which 
we find the velocity of the emission to be at 324± 15 km s _1 . 
Applying the if-correction, we find the projected veloc- 
ity of the centre of mass of the secondary star to be 
370 ± 40 km s . Combining this with existing parameters, 
we constrain the mass ratio to be 0.044±°'°°;| , and the mass 
function of the pulsar to be O^Iq'^ Mq. Hence, the mass 
of the pulsar is 0A8±°]{1/ sin 3 i Mq. 

Using various estimates of the binary inclination an- 
gle we find no evidence to suggest that the neutron star 
in SAX J1808.4— 3658 is more massive than the canonical 
value of 1.4 Mq, despite the fact that AMSPs are expected 
to have accreted a significant fraction of a solar mass in order 
to attain such rapid spin frequencies. 

An important next step in determining the mass of the 
neutron star in SAX J1808.4— 3658 is to confi rm the orbital 
inclin ation estimate from the Fe K line fitting l|Cackett et al.l 
2009). The outburst optical light curves appear to be domi- 
nated by the disc and not, for example, by an X-ray heated 
secondary that could be modelled to infer the inclination. 
In quiescence, however, a m odulation due to a heated sec- 
ondary is clearly visible (e.g. IDelove et alll2008l ). and more 
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detailed multicolour modelling here is probably the best 
route to confirming the inclination. 

Determining the masses of the neutron star in AMSPs 
is challenging for a number of reasons. However, because of 
the ability to accurately measure the orbital period and pro- 
jected velocity of the neutron star using the X-ray pulses, 
these AMSPs remain the most promising systems for deter- 
mining neutron star masses, helping to constrain the equa- 
tion of state which describes matter in these neturon stars. 
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